(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     81817,       2148]
NotebookOptionsPosition[     80124,       2086]
NotebookOutlinePosition[     80491,       2102]
CellTagsIndexPosition[     80448,       2099]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Setup MEH[k]", "Subsubsection",
 CellChangeTimes->{{3.570327517917389*^9, 3.57032751997117*^9}, {
  3.570327554725595*^9, 3.57032755584058*^9}, {3.570327787927311*^9, 
  3.5703277888672857`*^9}, {3.5703722949756393`*^9, 3.570372296323347*^9}, {
  3.592769661830781*^9, 3.592769666954217*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"c", " ", "=", " ", "299792.458"}], ";", " ", 
  RowBox[{"z", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h", "=", "0.70"}], ";", 
  RowBox[{"Omm", " ", "=", " ", "0.27"}], ";", 
  RowBox[{"Omb", "=", "0.04"}], ";", 
  RowBox[{"keq", " ", "=", " ", 
   FractionBox[
    RowBox[{"0.0746", "*", "Omm", "*", 
     SuperscriptBox["h", "2"]}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"2.728", "/", "2.7"}], ")"}], "2"]]}], ";", 
  RowBox[{"zdrag", "=", 
   RowBox[{"1291", "*", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Omm", "*", 
        RowBox[{"h", "^", "2"}]}], ")"}], "^", 
      RowBox[{"(", "0.251", ")"}]}], "/", 
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{"0.659", "*", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"Omm", "*", 
           RowBox[{"h", "^", "2"}]}], ")"}], "^", 
         RowBox[{"(", "0.828", ")"}]}]}]}], ")"}]}], "*", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"0.313", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"-", "0.419"}], ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"0.607", "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Omm", "*", 
             RowBox[{"h", "^", "2"}]}], ")"}], "^", "0.674"}]}]}], ")"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Omb", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"0.238", "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Omm", "*", 
             RowBox[{"h", "^", "2"}]}], ")"}], "^", "0.223"}]}], ")"}]}]}]}], 
     ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"zeq", "=", 
   RowBox[{"2.50", "*", 
    RowBox[{"10", "^", "4"}], "*", "Omm", "*", 
    RowBox[{
     RowBox[{"h", "^", "2"}], "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2.726", "/", "2.7"}], ")"}], "^", "4"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Rd", "=", 
    RowBox[{"31.5", "*", "Omb", "*", 
     RowBox[{"h", "^", "2"}], "*", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2.726", "/", "2.7"}], ")"}], "^", 
       RowBox[{"(", 
        RowBox[{"-", "4"}], ")"}]}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", "zdrag"}], ")"}], "/", 
        RowBox[{"10", "^", "3"}]}], ")"}]}]}]}], ";"}], 
  RowBox[{"(*", " ", 
   RowBox[{"R_drag", " ", "in", " ", 
    RowBox[{"tf_fit", ".", "c"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Req", "=", 
   RowBox[{"31.5", "*", "Omb", "*", 
    RowBox[{"h", "^", "2"}], "*", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2.726", "/", "2.7"}], ")"}], "^", 
      RowBox[{"(", 
       RowBox[{"-", "4"}], ")"}]}], "/", 
     RowBox[{"(", 
      RowBox[{"zeq", "/", 
       RowBox[{"10", "^", "3"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 

 RowBox[{
  RowBox[{
   RowBox[{"kSilk", "=", 
    RowBox[{"1.6", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Omb", "*", 
        RowBox[{"h", "^", "2"}]}], ")"}], "^", 
      RowBox[{"(", "0.52", ")"}]}], "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Omm", "*", 
        RowBox[{"h", "^", "2"}]}], ")"}], "^", 
      RowBox[{"(", "0.73", ")"}]}], "*", 
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"10.4", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"-", "0.95"}], ")"}]}]}], ")"}]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"in", " ", 
    RowBox[{"Mpc", "^", 
     RowBox[{"(", 
      RowBox[{"-", "1"}], ")"}]}]}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"xx", "=", 
    RowBox[{
     RowBox[{
      RowBox[{"2", "/", "3"}], "/", "keq"}], "*", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"6", "/", "Req"}], "]"}], "*", 
     RowBox[{"Log", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Sqrt", "[", 
          RowBox[{"1", "+", "Rd"}], "]"}], "+", 
         RowBox[{"Sqrt", "[", 
          RowBox[{"Rd", "+", "Req"}], "]"}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"Sqrt", "[", "Req", "]"}]}], ")"}]}], "]"}]}]}], ";"}], 
  StyleBox[
   RowBox[{"(*", " ", "sound_horizon", " ", "*)"}],
   FontColor->GrayLevel[0]]}], "\n", 
 RowBox[{
  RowBox[{"\[Beta]c", "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{
       RowBox[{"0.944", "/", 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"458", "*", "Omm", "*", 
             RowBox[{"h", "^", "2"}]}], ")"}], "^", 
           RowBox[{"(", 
            RowBox[{"-", "0.708"}], ")"}]}]}], ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"Omm", "-", "Omb"}], ")"}], "/", "Omm"}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.395", "*", "Omm", "*", 
              RowBox[{"h", "^", "2"}]}], ")"}], "^", 
            RowBox[{"(", 
             RowBox[{"-", "0.0266"}], ")"}]}], ")"}]}], "-", "1"}], ")"}]}]}],
      ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Beta]b", "=", 
   RowBox[{"0.5", "+", 
    RowBox[{"Omb", "/", "Omm"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"3", "-", 
       RowBox[{"2", "*", 
        RowBox[{"Omb", "/", "Omm"}]}]}], ")"}], "*", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"17.2", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", "2"}], "+", "1"}], 
      "]"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Alpha]c", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"46.9", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", "0.670", ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"32.1", "*", "Omm", "*", 
            RowBox[{"h", "^", "2"}]}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"-", "0.532"}], ")"}]}]}], ")"}]}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "Omb"}], "/", "Omm"}], ")"}]}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"12.0", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", "0.424", ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"45.0", "*", "Omm", "*", 
            RowBox[{"h", "^", "2"}]}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"-", "0.582"}], ")"}]}]}], ")"}]}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{"-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Omb", "/", "Omm"}], ")"}], "^", "3"}]}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Alpha]b", "=", 
   RowBox[{"2.07", "*", "keq", "*", "xx", "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", "Rd"}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "3"}], "/", "4"}], ")"}]}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", "zeq"}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{"1", "+", "zdrag"}], ")"}]}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "6"}], "*", 
       RowBox[{"Sqrt", "[", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "+", "zeq"}], ")"}], "/", 
          RowBox[{"(", 
           RowBox[{"1", "+", "zdrag"}], ")"}]}]}], "]"}]}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "+", 
         RowBox[{"3", "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "+", "zeq"}], ")"}], "/", 
           RowBox[{"(", 
            RowBox[{"1", "+", "zdrag"}], ")"}]}]}]}], ")"}], "*", 
       RowBox[{"Log", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Sqrt", "[", 
            RowBox[{"1", "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1", "+", "zeq"}], ")"}], "/", 
              RowBox[{"(", 
               RowBox[{"1", "+", "zdrag"}], ")"}]}]}], "]"}], "+", "1"}], 
          ")"}], "/", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Sqrt", "[", 
            RowBox[{"1", "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1", "+", "zeq"}], ")"}], "/", 
              RowBox[{"(", 
               RowBox[{"1", "+", "zdrag"}], ")"}]}]}], "]"}], "-", "1"}], 
          ")"}]}], "]"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TransferEH", "[", "k_", "]"}], ":=", " ", 
   RowBox[{"Piecewise", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Omm", "-", "Omb"}], ")"}], "/", "Omm"}], "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "/", 
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"k", "*", 
                    RowBox[{"xx", "/", "5.4"}]}], ")"}], "^", "4"}]}], 
                ")"}]}], ")"}], "*", 
             RowBox[{
              RowBox[{"Log", "[", 
               RowBox[{"E", "+", 
                RowBox[{"1.8", "*", "\[Beta]c", "*", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
               "]"}], "/", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Log", "[", 
                 RowBox[{"E", "+", 
                  RowBox[{"1.8", "*", "\[Beta]c", "*", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                 "]"}], "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"14.2", "+", 
                   RowBox[{"386", "/", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    RowBox[{"69.9", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                    "1.08"}]}]}], ")"}]}]}], ")"}], "*", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                  "2"}]}]}], ")"}]}]}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "-", 
               RowBox[{"(", 
                RowBox[{"1", "/", 
                 RowBox[{"(", 
                  RowBox[{"1", "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"k", "*", 
                    RowBox[{"xx", "/", "5.4"}]}], ")"}], "^", "4"}]}], 
                  ")"}]}], ")"}]}], ")"}], "*", 
             RowBox[{
              RowBox[{"Log", "[", 
               RowBox[{"E", "+", 
                RowBox[{"1.8", "*", "\[Beta]c", "*", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
               "]"}], "/", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Log", "[", 
                 RowBox[{"E", "+", 
                  RowBox[{"1.8", "*", "\[Beta]c", "*", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                 "]"}], "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"14.2", "/", "\[Alpha]c"}], "+", 
                   RowBox[{"386", "/", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    RowBox[{"69.9", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                    "1.08"}]}]}], ")"}]}]}], ")"}], "*", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                  "2"}]}]}], ")"}]}]}]}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"Omb", "/", "Omm"}], "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Log", "[", 
                RowBox[{"E", "+", 
                 RowBox[{"1.8", "*", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                "]"}], "/", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"Log", "[", 
                  RowBox[{"E", "+", 
                   RowBox[{"1.8", "*", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                  "]"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"14.2", "+", 
                    RowBox[{"386", "/", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    RowBox[{"69.9", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                    RowBox[{"(", "1.08", ")"}]}]}]}], ")"}]}]}], ")"}], "*", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                   "2"}]}]}], ")"}]}], ")"}], "/", 
             RowBox[{"(", 
              RowBox[{"1", "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"k", "*", 
                  RowBox[{"xx", "/", "5.2"}]}], ")"}], "^", "2"}]}], ")"}]}], 
            "+", 
            RowBox[{"\[Alpha]b", "*", 
             RowBox[{
              RowBox[{"Exp", "[", 
               RowBox[{"-", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"k", "/", "kSilk"}], ")"}], "^", "1.4"}]}], "]"}], 
              "/", 
              RowBox[{"(", 
               RowBox[{"1", "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"\[Beta]b", "/", "k"}], "/", "xx"}], ")"}], "^", 
                 "3"}]}], ")"}]}]}]}], ")"}], 
          RowBox[{
           RowBox[{"Sin", "[", 
            RowBox[{"k", "*", 
             RowBox[{"xx", "/", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"8.41", "*", 
                    RowBox[{
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"Omm", "*", 
                    RowBox[{"h", "^", "2"}]}], ")"}], "^", 
                    RowBox[{"(", "0.435", ")"}]}], "/", "k"}], "/", "xx"}]}], 
                   ")"}], "^", "3"}]}], ")"}], "^", 
               RowBox[{"(", 
                RowBox[{"1", "/", "3"}], ")"}]}]}]}], "]"}], "/", 
           RowBox[{"(", 
            RowBox[{"k", "*", 
             RowBox[{"xx", "/", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"8.41", "*", 
                    RowBox[{
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"Omm", "*", 
                    RowBox[{"h", "^", "2"}]}], ")"}], "^", 
                    RowBox[{"(", "0.435", ")"}]}], "/", "k"}], "/", "xx"}]}], 
                   ")"}], "^", "3"}]}], ")"}], "^", 
               RowBox[{"(", 
                RowBox[{"1", "/", "3"}], ")"}]}]}]}], ")"}]}]}]}], ",", 
        RowBox[{"k", ">", "0"}]}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"1.0", ",", " ", 
        RowBox[{"k", "\[Equal]", "0"}]}], "}"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Hubble", "[", "zz_", "]"}], ":=", 
    RowBox[{"100", " ", "h", " ", 
     SqrtBox[
      RowBox[{
       RowBox[{"Omm", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"1", "+", "zz"}], ")"}], "3"]}], "+", 
       RowBox[{"(", 
        RowBox[{"1", "-", "Omm"}], ")"}]}]]}]}], ";"}], "\n", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"don", "'"}], "t", " ", "forget", " ", "below", " ", "is", " ", 
    "only", " ", "valid", " ", "for", " ", "a", " ", "flat", " ", 
    "\[CapitalLambda]CDM", " ", "or", " ", "CDM", " ", "universe"}], 
   "*)"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Growth", "[", "zz_", "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{"5", " ", "Omm"}], "2"], " ", 
    FractionBox[
     RowBox[{"Hubble", "[", "zz", "]"}], 
     RowBox[{"100", " ", "h"}]], 
    RowBox[{"NIntegrate", "[", 
     RowBox[{
      FractionBox[
       RowBox[{"(", 
        RowBox[{"1", "+", "y"}], ")"}], 
       SuperscriptBox[
        RowBox[{"(", 
         FractionBox[
          RowBox[{"Hubble", "[", "y", "]"}], 
          RowBox[{"100", " ", "h"}]], ")"}], "3"]], ",", 
      RowBox[{"{", 
       RowBox[{"y", ",", "zz", ",", "\[Infinity]"}], "}"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"MEH", "[", "k_", "]"}], ":=", 
    RowBox[{
     RowBox[{
      RowBox[{"2.", "/", "3."}], "/", "Omm"}], "*", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"c", "/", 
        RowBox[{"(", 
         RowBox[{"100", " ", "h"}], ")"}]}], ")"}], "2"], "*", 
     RowBox[{"TransferEH", "[", "k", "]"}], "*", 
     SuperscriptBox["k", "2"], "*", 
     RowBox[{"Growth", "[", "z", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"Simple", " ", 
    RowBox[{"Transfer", "[", "k", "]"}]}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Omb", " ", "=", " ", "0.04"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"baryon", " ", "density"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "=", 
   RowBox[{"Omm", "*", "h", "*", 
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", "Omb"}], 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{
         RowBox[{"Sqrt", "[", 
          RowBox[{"2", "*", "h"}], "]"}], "/", "Omm"}]}], ")"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Transfer", "[", "k_", "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"Log", "[", 
      RowBox[{"1", "+", 
       RowBox[{"2.34", 
        RowBox[{"(", 
         RowBox[{"k", "/", 
          RowBox[{"(", 
           RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}]}], "]"}], 
     "/", 
     RowBox[{"(", 
      RowBox[{"2.34", "*", 
       RowBox[{"(", 
        RowBox[{"k", "/", 
         RowBox[{"(", 
          RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}]}], 
    "*", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{"3.89", "*", 
        RowBox[{"(", 
         RowBox[{"k", "/", 
          RowBox[{"(", 
           RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"16.1", "*", 
          RowBox[{"(", 
           RowBox[{"k", "/", 
            RowBox[{"(", 
             RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}], 
        "2"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"5.46", "*", 
          RowBox[{"(", 
           RowBox[{"k", "/", 
            RowBox[{"(", 
             RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}], 
        "3"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"6.71", "*", 
          RowBox[{"(", 
           RowBox[{"k", "/", 
            RowBox[{"(", 
             RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}], 
        "4"]}], ")"}], 
     RowBox[{"-", "0.25"}]]}]}], ";"}]}], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.563542383318079*^9, 3.563542460983938*^9}, 
   3.563542493723468*^9, {3.563542586887805*^9, 3.563542628749124*^9}, {
   3.563543087042822*^9, 3.563543090748076*^9}, {3.563543374436173*^9, 
   3.563543422198125*^9}, {3.5644021704266863`*^9, 3.564402172967256*^9}, {
   3.564405767539486*^9, 3.564405768052952*^9}, {3.564413235949687*^9, 
   3.564413236853958*^9}, {3.564682427610278*^9, 3.564682436448017*^9}, {
   3.5646824745485983`*^9, 3.564682480322901*^9}, {3.5648343000988703`*^9, 
   3.56483430413787*^9}, {3.564850105889791*^9, 3.56485010797779*^9}, {
   3.566426830196505*^9, 3.566426831159525*^9}, 3.570327537517303*^9, {
   3.570327712775303*^9, 3.570327714117038*^9}, {3.5703768187944307`*^9, 
   3.570376821241983*^9}, {3.5715755872881193`*^9, 3.571575636954029*^9}, {
   3.5885027348150806`*^9, 3.588502735071302*^9}, {3.588502813821426*^9, 
   3.588502814028502*^9}, {3.5885032878929996`*^9, 3.588503339433796*^9}, {
   3.588503406456326*^9, 3.588503406639093*^9}, {3.588503465987542*^9, 
   3.588503471923272*^9}, {3.5922255614844513`*^9, 3.592225562129325*^9}, {
   3.5927643164655037`*^9, 3.592764317763259*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MEH", "[", 
  RowBox[{"0.01", "*", "h"}], "]"}]], "Input",
 CellChangeTimes->{{3.592152851087934*^9, 3.592152852570692*^9}, {
  3.592152889472789*^9, 3.592152895286504*^9}}],

Cell[BoxData["1275.1632126042107`"], "Output",
 CellChangeTimes->{3.592152895936891*^9, 3.592223409910748*^9, 
  3.592772110419016*^9, 3.592772234865951*^9, 3.5929147545551453`*^9, 
  3.5931914766945868`*^9, 3.593385026671564*^9, 3.593387544136036*^9, 
  3.593443816729603*^9, 3.593443973908484*^9, 3.593445028950808*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LogLogPlot", "[", 
  RowBox[{
   RowBox[{
    SuperscriptBox[
     RowBox[{"MEH", "[", 
      RowBox[{"k", "*", "h"}], "]"}], "2"], "*", 
    RowBox[{"Pvar", "[", "k", "]"}]}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"k", ",", " ", "0.001", ",", "1"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.592223642355777*^9, 3.592223770239456*^9}, {
   3.5922238479599323`*^9, 3.5922238482505207`*^9}, {3.592223900806583*^9, 
   3.592223902446669*^9}, 3.592223973952518*^9}],

Cell[BoxData[
 GraphicsBox[{},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{-6.907755278982137, -0.6931471805599453},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Exp[
       Part[#, 1]], 
      Exp[
       Part[#, 2]]}& ), "CopiedValueFunction" -> ({
      Exp[
       Part[#, 1]], 
      Exp[
       Part[#, 2]]}& )},
  FrameTicks->{{{{0., 
       FormBox[
        TagBox[
         InterpretationBox["\"1.0\"", 1., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {2.302585092994046, 
       FormBox[
        TagBox[
         InterpretationBox["\"10.0\"", 10., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.6094379124341003`, 
       FormBox[
        TagBox[
         InterpretationBox["\"5.0\"", 5., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {3.912023005428146, 
       FormBox[
        TagBox[
         InterpretationBox["\"50.0\"", 50., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0.6931471805599453, 
       FormBox[
        TagBox[
         InterpretationBox["\"2.0\"", 2., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {2.995732273553991, 
       FormBox[
        TagBox[
         InterpretationBox["\"20.0\"", 20., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.0986122886681098`, 
       FormBox[
        TagBox[
         InterpretationBox["\"3.0\"", 3., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {3.4011973816621555`, 
       FormBox[
        TagBox[
         InterpretationBox["\"30.0\"", 30., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0.4054651081081644, 
       FormBox[
        TagBox[
         InterpretationBox["\"1.5\"", 1.5, AutoDelete -> True], NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {2.70805020110221, 
       FormBox[
        TagBox[
         InterpretationBox["\"15.0\"", 15., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.9459101490553132`, 
       FormBox[
        TagBox[
         InterpretationBox["\"7.0\"", 7., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {4.248495242049359, 
       FormBox[
        TagBox[
         InterpretationBox["\"70.0\"", 70., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0.09531017980432493, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.1823215567939548, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.26236426446749106`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.336472236621213, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.47000362924573563`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5306282510621705, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5877866649021191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6418538861723948, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.791759469228055, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.0794415416798357`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.1972245773362196`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.3978952727983707`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.4849066497880004`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.5649493574615367`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.6390573296152584`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.772588722239781, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.833213344056216, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.8903717578961645`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.9444389791664403`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {3.6888794541139363`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {4.0943445622221, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}, {{0., 
       FormBox["\"\"", TraditionalForm]}, {2.302585092994046, 
       FormBox["\"\"", TraditionalForm]}, {1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm]}, {3.912023005428146, 
       FormBox["\"\"", TraditionalForm]}, {0.6931471805599453, 
       FormBox["\"\"", TraditionalForm]}, {2.995732273553991, 
       FormBox["\"\"", TraditionalForm]}, {1.0986122886681098`, 
       FormBox["\"\"", TraditionalForm]}, {3.4011973816621555`, 
       FormBox["\"\"", TraditionalForm]}, {0.4054651081081644, 
       FormBox["\"\"", TraditionalForm]}, {2.70805020110221, 
       FormBox["\"\"", TraditionalForm]}, {1.9459101490553132`, 
       FormBox["\"\"", TraditionalForm]}, {4.248495242049359, 
       FormBox["\"\"", TraditionalForm]}, {0.09531017980432493, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.1823215567939548, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.26236426446749106`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.336472236621213, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.47000362924573563`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5306282510621705, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5877866649021191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6418538861723948, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.791759469228055, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.0794415416798357`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.1972245773362196`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.3978952727983707`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.4849066497880004`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.5649493574615367`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.6390573296152584`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.772588722239781, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.833213344056216, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.8903717578961645`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.9444389791664403`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {3.6888794541139363`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {4.0943445622221, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}}, {{{-6.907755278982137, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.001\"", 0.001, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-5.298317366548036, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.005\"", 0.005, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-4.605170185988092, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.010\"", 0.01, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-2.995732273553991, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.050\"", 0.05, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-2.302585092994046, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.100\"", 0.1, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-0.6931471805599453, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.500\"", 0.5, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0., 
       FormBox[
        TagBox[
         InterpretationBox["\"1.000\"", 1., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.6094379124341003`, 
       FormBox[
        TagBox[
         InterpretationBox["\"5.000\"", 5., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-6.214608098422191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.809142990314028, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.521460917862246, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.115995809754082, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.961845129926823, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.8283137373023015`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.710530701645918, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.912023005428146, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.506557897319982, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.2188758248682006`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.8134107167600364`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.6592600369327783`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.5257286443082556`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.4079456086518722`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.203972804325936, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.916290731874155, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.5108256237659907, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.35667494393873245`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.2231435513142097, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.10536051565782628`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6931471805599453, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.0986122886681098`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}, {{-6.907755278982137, 
       FormBox["\"\"", TraditionalForm]}, {-5.298317366548036, 
       FormBox["\"\"", TraditionalForm]}, {-4.605170185988092, 
       FormBox["\"\"", TraditionalForm]}, {-2.995732273553991, 
       FormBox["\"\"", TraditionalForm]}, {-2.302585092994046, 
       FormBox["\"\"", TraditionalForm]}, {-0.6931471805599453, 
       FormBox["\"\"", TraditionalForm]}, {0., 
       FormBox["\"\"", TraditionalForm]}, {1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm]}, {-6.214608098422191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.809142990314028, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.521460917862246, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.115995809754082, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.961845129926823, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.8283137373023015`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.710530701645918, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.912023005428146, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.506557897319982, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.2188758248682006`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.8134107167600364`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.6592600369327783`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.5257286443082556`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.4079456086518722`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.203972804325936, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.916290731874155, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.5108256237659907, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.35667494393873245`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.2231435513142097, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.10536051565782628`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6931471805599453, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.0986122886681098`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}}},
  GridLines->{None, None},
  Method->{},
  PlotRange->{{-6.907755278982137, 0}, {-0.6931471805599453, 
    0.6931471805599453}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]},
  Ticks->{{{-6.907755278982137, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.001\"", 0.001, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-5.298317366548036, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.005\"", 0.005, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-4.605170185988092, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.010\"", 0.01, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-2.995732273553991, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.050\"", 0.05, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-2.302585092994046, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.100\"", 0.1, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-0.6931471805599453, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.500\"", 0.5, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0., 
      FormBox[
       TagBox[
        InterpretationBox["\"1.000\"", 1., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.6094379124341003`, 
      FormBox[
       TagBox[
        InterpretationBox["\"5.000\"", 5., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-6.214608098422191, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-5.809142990314028, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-5.521460917862246, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-5.115995809754082, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-4.961845129926823, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-4.8283137373023015`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-4.710530701645918, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-3.912023005428146, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-3.506557897319982, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-3.2188758248682006`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.8134107167600364`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.6592600369327783`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.5257286443082556`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.4079456086518722`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-1.6094379124341003`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-1.203972804325936, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.916290731874155, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.5108256237659907, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.35667494393873245`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.2231435513142097, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.10536051565782628`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.6931471805599453, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.0986122886681098`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.3862943611198906`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}}, {{0., 
      FormBox[
       TagBox[
        InterpretationBox["\"1.0\"", 1., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {2.302585092994046, 
      FormBox[
       TagBox[
        InterpretationBox["\"10.0\"", 10., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.6094379124341003`, 
      FormBox[
       TagBox[
        InterpretationBox["\"5.0\"", 5., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {3.912023005428146, 
      FormBox[
       TagBox[
        InterpretationBox["\"50.0\"", 50., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0.6931471805599453, 
      FormBox[
       TagBox[
        InterpretationBox["\"2.0\"", 2., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {2.995732273553991, 
      FormBox[
       TagBox[
        InterpretationBox["\"20.0\"", 20., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.0986122886681098`, 
      FormBox[
       TagBox[
        InterpretationBox["\"3.0\"", 3., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {3.4011973816621555`, 
      FormBox[
       TagBox[
        InterpretationBox["\"30.0\"", 30., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0.4054651081081644, 
      FormBox[
       TagBox[
        InterpretationBox["\"1.5\"", 1.5, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {2.70805020110221, 
      FormBox[
       TagBox[
        InterpretationBox["\"15.0\"", 15., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.9459101490553132`, 
      FormBox[
       TagBox[
        InterpretationBox["\"7.0\"", 7., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {4.248495242049359, 
      FormBox[
       TagBox[
        InterpretationBox["\"70.0\"", 70., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0.09531017980432493, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.1823215567939548, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.26236426446749106`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.336472236621213, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.47000362924573563`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.5306282510621705, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.5877866649021191, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.6418538861723948, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.3862943611198906`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.791759469228055, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.0794415416798357`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.1972245773362196`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.3978952727983707`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.4849066497880004`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.5649493574615367`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.6390573296152584`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.772588722239781, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.833213344056216, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.8903717578961645`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.9444389791664403`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {3.6888794541139363`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {4.0943445622221, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}}}]], "Output",
 CellChangeTimes->{{3.592223720814096*^9, 3.5922237762299767`*^9}, 
   3.5922238545593777`*^9, 3.5922239074637327`*^9, 3.592223979082264*^9, 
   3.592772113604527*^9, 3.5927722376036997`*^9, 3.5929147573631077`*^9, 
   3.5931914796729517`*^9, 3.593385029277939*^9, 3.593387546097588*^9, 
   3.59344382090629*^9, 3.593443979026782*^9, 3.593445032682391*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"D", "[", 
  RowBox[{
   RowBox[{"TopHat", "[", 
    RowBox[{"k", ",", " ", "RR"}], "]"}], ",", " ", "RR"}], "]"}]], "Input",
 CellChangeTimes->{{3.592869572639958*^9, 3.5928695813045273`*^9}}],

Cell[BoxData[
 RowBox[{
  SuperscriptBox["TopHat", 
   TagBox[
    RowBox[{"(", 
     RowBox[{"0", ",", "1"}], ")"}],
    Derivative],
   MultilineFunction->None], "[", 
  RowBox[{"k", ",", "RR"}], "]"}]], "Output",
 CellChangeTimes->{3.592869581704977*^9, 3.592914757432687*^9, 
  3.5931914797866697`*^9, 3.5933850293324757`*^9, 3.5933875461475677`*^9, 
  3.593443821359991*^9, 3.5934439797877502`*^9, 3.593445033132846*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Setup Integrals", "Subsection",
 CellChangeTimes->{{3.584276308436615*^9, 3.584276368172073*^9}, {
   3.58427647538827*^9, 3.584276476437305*^9}, 3.584276521317811*^9, {
   3.5885033280710697`*^9, 3.5885033316502867`*^9}, {3.59206949221636*^9, 
   3.5920694960873213`*^9}, {3.592769687852216*^9, 3.5927696906303167`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"k0", "=", 
    RowBox[{"0.002", "/", "h"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"In", " ", 
    RowBox[{"h", "/", "Mpc"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ns", "=", "0.96"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Lx", "=", "2400"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"In", " ", 
    RowBox[{"Mpc", "/", "h"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Np", "=", 
   SuperscriptBox["1024", "3"]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Mp", " ", "=", " ", 
    RowBox[{"9.648", " ", 
     SuperscriptBox["10", "11"]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"mass", " ", "of", " ", "a", " ", "particle", " ", "in", " ", 
    SuperscriptBox["h", 
     RowBox[{"-", "1"}]], 
    SubscriptBox["M", "sun"]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Rho]bar", " ", "=", " ", 
    FractionBox[
     RowBox[{"Np", "*", "Mp"}], 
     SuperscriptBox[
      RowBox[{"(", "Lx", ")"}], "3"]]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"in", " ", 
    SuperscriptBox["h", "2"], 
    RowBox[{
     SubscriptBox["M", "sun"], "/", 
     SuperscriptBox["Mpc", "3"]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Mass", "[", "R_", "]"}], ":=", " ", 
    RowBox[{
     FractionBox["4", "3"], "\[Pi]", " ", 
     SuperscriptBox["R", "3"], " ", "*", "\[Rho]bar"}]}], ";"}], "   ", 
  RowBox[{"(*", " ", 
   RowBox[{"in", " ", 
    SuperscriptBox["h", 
     RowBox[{"-", "1"}]], 
    SubscriptBox["M", "sun"]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Aphi", "=", 
   RowBox[{"8.23", " ", 
    SuperscriptBox["10", 
     RowBox[{"-", "10"}]]}]}], ";", 
  RowBox[{"Aphi", "=", "1."}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"kmin", "=", 
    RowBox[{"2", "*", 
     RowBox[{"\[Pi]", "/", "Lx"}]}]}], ";", " ", 
   RowBox[{"kmax", "=", 
    RowBox[{"1024", "*", "kmin"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"in", " ", 
    RowBox[{"h", "/", "Mpc"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Pvar", "[", "k_", "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"k", ">=", "kmin"}], ",", " ", 
     RowBox[{"Aphi", "*", 
      RowBox[{
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"k", "/", "k0"}], ")"}], 
        RowBox[{"ns", "-", "1"}]], "/", 
       SuperscriptBox["k", "3"]}]}], ",", " ", "0."}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TopHat", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
   FractionBox[
    RowBox[{
     RowBox[{"3", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"k", " ", "R"}], "]"}]}], "-", 
     RowBox[{"3", 
      RowBox[{"(", 
       RowBox[{"k", " ", "R"}], ")"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"k", " ", "R"}], "]"}]}]}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"k", " ", "R"}], ")"}], "3"]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TopHatD", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{"3", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"k", " ", "R"}], "]"}]}], 
     RowBox[{"k", " ", 
      SuperscriptBox["R", "2"]}]], "-", 
    FractionBox[
     RowBox[{"3", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "3"}], " ", "k", " ", "R", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"k", " ", "R"}], "]"}]}], "+", 
        RowBox[{"3", " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"k", " ", "R"}], "]"}]}]}], ")"}]}], 
     RowBox[{
      SuperscriptBox["k", "3"], " ", 
      SuperscriptBox["R", "4"]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Alpha]", "[", "k_", "]"}], ":=", 
   RowBox[{"MEH", "[", 
    RowBox[{"k", "*", "h"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"\[Alpha]M", "[", 
     RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"MEH", "[", 
      RowBox[{"k", "*", "h"}], "]"}], "*", 
     RowBox[{"TopHat", "[", 
      RowBox[{"k", ",", "R"}], "]"}]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"MEH", "[", "k", "]"}], " ", "takes", " ", "k", " ", "in", " ", 
    SuperscriptBox["Mpc", 
     RowBox[{"-", "1"}]]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"CosTheta", "[", 
    RowBox[{"\[Phi]1_", ",", "\[Phi]2_", ",", "\[Mu]1_", ",", "\[Mu]2_"}], 
    "]"}], ":=", 
   RowBox[{
    RowBox[{"\[Mu]1", "*", "\[Mu]2"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"Cos", "[", "\[Phi]1", "]"}], "*", 
        RowBox[{"Cos", "[", "\[Phi]2", "]"}]}], "+", 
       RowBox[{
        RowBox[{"Sin", "[", "\[Phi]1", "]"}], "*", 
        RowBox[{"Sin", "[", "\[Phi]2", "]"}]}]}], ")"}], "*", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"\[Mu]1", "^", "2"}]}], ")"}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"\[Mu]2", "^", "2"}]}], ")"}]}], "]"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I32Value", "[", 
    RowBox[{
    "k_", ",", "R_", ",", " ", "q2_", ",", " ", "\[Mu]2_", ",", " ", 
     "\[Phi]2_", ",", " ", "q1_", ",", " ", "\[Mu]1_", ",", " ", "\[Phi]1_"}],
     "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["q2", "2"], "*", 
       RowBox[{"\[Alpha]M", "[", 
        RowBox[{"q2", ",", "R"}], "]"}], "*", 
       RowBox[{"\[Alpha]M", "[", 
        RowBox[{
         RowBox[{"Sqrt", "[", 
          RowBox[{
           SuperscriptBox["q2", "2"], "+", 
           SuperscriptBox["k", "2"], "+", 
           RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], ",", 
         "R"}], "]"}], "*", 
       SuperscriptBox["q1", "2"], "*", 
       RowBox[{"Pvar", "[", "q1", "]"}], 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q2", "2"], "+", 
          SuperscriptBox["q1", "2"], "-", 
          RowBox[{"2", "*", "q2", "*", "q1", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{
            "\[Phi]1", ",", " ", "\[Phi]2", ",", " ", "\[Mu]1", ",", " ", 
             "\[Mu]2"}], "]"}]}]}], "]"}], "]"}], "*", 
       RowBox[{"(", 
        RowBox[{"Pvar", "[", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           SuperscriptBox["k", "2"], "+", 
           SuperscriptBox["q1", "2"], "+", 
           RowBox[{"2", "*", "k", "*", "q1", "*", "\[Mu]1"}]}], "]"}], "]"}], 
        ")"}]}], ")"}], "/", "8."}], "/", 
    SuperscriptBox["\[Pi]", "2"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I31Value", "[", 
    RowBox[{"k_", ",", "R_", ",", "q1_", ",", "\[Mu]1_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox["q1", "2"], "*", 
    RowBox[{"\[Alpha]M", "[", 
     RowBox[{"q1", ",", "R"}], "]"}], "*", 
    RowBox[{"\[Alpha]M", "[", 
     RowBox[{
      RowBox[{"Sqrt", "[", 
       RowBox[{
        SuperscriptBox["q1", "2"], "+", 
        SuperscriptBox["k", "2"], "+", 
        RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], ",", 
      "R"}], "]"}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Pvar", "[", "q1", "]"}], "*", 
         RowBox[{"Pvar", "[", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], 
          "]"}]}], ")"}], "+", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Pvar", "[", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], "]"}],
          "*", 
         RowBox[{"Pvar", "[", "k", "]"}]}], ")"}], "+", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Pvar", "[", "k", "]"}], "*", 
         RowBox[{"Pvar", "[", "q1", "]"}]}], ")"}]}], ")"}], "/", "4.0"}]}]}],
   ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"ISigSq", "[", 
     RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"\[Alpha]M", "[", 
       RowBox[{"k", ",", "R"}], "]"}], "2"], "*", 
     RowBox[{"Pvar", "[", "k", "]"}], "*", "k", "*", "k"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "for", " ", "SigmaMSq", " ", "computation", " ", "to", " ", "get", " ", 
     "A", " ", "such", " ", "that", " ", 
     SubscriptBox["\[Sigma]", "8"]}], "=", "0.8"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{",", " ", 
      RowBox[{"{", 
       RowBox[{"q2", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Mu]2", ",", " ", 
        RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Phi]2", ",", " ", "0", ",", " ", 
        RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"q1", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Mu]1", ",", " ", 
        RowBox[{"-", "1"}], ",", " ", "1"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Phi]1", ",", " ", "0", ",", " ", 
        RowBox[{"2", "*", "\[Pi]"}]}], "}"}]}], "]"}], ";"}], "*)"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"I31DLog", "[", 
     RowBox[{"k_", ",", "R_", ",", "q1_", ",", "\[Mu]1_"}], "]"}], ":=", 
    RowBox[{
     SuperscriptBox["q1", "2"], "*", 
     RowBox[{"\[Alpha]", "[", "q1", "]"}], "*", 
     RowBox[{"\[Alpha]", "[", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        SuperscriptBox["q1", "2"], "+", 
        SuperscriptBox["k", "2"], "+", 
        RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], "]"}], 
     "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"TopHatD", "[", 
         RowBox[{"q1", ",", "R"}], "]"}], "*", 
        RowBox[{"TopHat", "[", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], ",", 
          "R"}], "]"}]}], "+", 
       RowBox[{
        RowBox[{"TopHat", "[", 
         RowBox[{"q1", ",", "R"}], "]"}], "*", 
        RowBox[{"TopHatD", "[", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], ",", 
          "R"}], "]"}]}]}], ")"}], "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Pvar", "[", "q1", "]"}], "*", 
          RowBox[{"Pvar", "[", 
           RowBox[{"Sqrt", "[", 
            RowBox[{
             SuperscriptBox["q1", "2"], "+", 
             SuperscriptBox["k", "2"], "+", 
             RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], 
           "]"}]}], ")"}], "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Pvar", "[", 
           RowBox[{"Sqrt", "[", 
            RowBox[{
             SuperscriptBox["q1", "2"], "+", 
             SuperscriptBox["k", "2"], "+", 
             RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], 
           "]"}], "*", 
          RowBox[{"Pvar", "[", "k", "]"}]}], ")"}], "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Pvar", "[", "k", "]"}], "*", 
          RowBox[{"Pvar", "[", "q1", "]"}]}], ")"}]}], ")"}], "/", 
      "4.0"}]}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I32DLog", "[", 
    RowBox[{
    "k_", ",", "R_", ",", " ", "q2_", ",", " ", "\[Mu]2_", ",", " ", 
     "\[Phi]2_", ",", " ", "q1_", ",", " ", "\[Mu]1_", ",", " ", "\[Phi]1_"}],
     "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["q2", "2"], "*", 
       RowBox[{"\[Alpha]", "[", "q2", "]"}], "*", 
       RowBox[{"\[Alpha]", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q2", "2"], "+", 
          SuperscriptBox["k", "2"], "+", 
          RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], "]"}], 
       "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"TopHatD", "[", 
           RowBox[{"q2", ",", "R"}], "]"}], "*", 
          RowBox[{"TopHat", "[", 
           RowBox[{
            RowBox[{"Sqrt", "[", 
             RowBox[{
              SuperscriptBox["q2", "2"], "+", 
              SuperscriptBox["k", "2"], "+", 
              RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], ",",
             "R"}], "]"}]}], "+", 
         RowBox[{
          RowBox[{"TopHat", "[", 
           RowBox[{"q2", ",", "R"}], "]"}], "*", 
          RowBox[{"TopHatD", "[", 
           RowBox[{
            RowBox[{"Sqrt", "[", 
             RowBox[{
              SuperscriptBox["q2", "2"], "+", 
              SuperscriptBox["k", "2"], "+", 
              RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], ",",
             "R"}], "]"}]}]}], ")"}], "*", 
       SuperscriptBox["q1", "2"], "*", 
       RowBox[{"Pvar", "[", "q1", "]"}], 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q2", "2"], "+", 
          SuperscriptBox["q1", "2"], "-", 
          RowBox[{"2", "*", "q2", "*", "q1", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{
            "\[Phi]1", ",", " ", "\[Phi]2", ",", " ", "\[Mu]1", ",", " ", 
             "\[Mu]2"}], "]"}]}]}], "]"}], "]"}], "*", 
       RowBox[{"(", 
        RowBox[{"Pvar", "[", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           SuperscriptBox["k", "2"], "+", 
           SuperscriptBox["q1", "2"], "+", 
           RowBox[{"2", "*", "k", "*", "q1", "*", "\[Mu]1"}]}], "]"}], "]"}], 
        ")"}]}], ")"}], "/", "8."}], "/", 
    SuperscriptBox["\[Pi]", "2"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I42Value", "[", 
    RowBox[{
    "k_", ",", "R_", ",", " ", "q2_", ",", " ", "\[Mu]2_", ",", " ", 
     "\[Phi]2_", ",", " ", "q1_", ",", " ", "\[Mu]1_", ",", " ", "\[Phi]1_", 
     ",", " ", "p_", ",", " ", "\[Mu]_", ",", " ", "\[Phi]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["q1", "2"], "*", 
       SuperscriptBox["q2", "2"], "*", 
       SuperscriptBox["p", "2"], "*", 
       SuperscriptBox[
        RowBox[{"\[Alpha]M", "[", 
         RowBox[{"q1", ",", "R"}], "]"}], "2"], "*", 
       SuperscriptBox[
        RowBox[{"\[Alpha]M", "[", 
         RowBox[{"q2", ",", "R"}], "]"}], "2"], "*", 
       RowBox[{"Pvar", "[", "p", "]"}], "*", 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q1", "2"], "-", 
          RowBox[{"2", "*", "q1", "*", "p", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]1", ",", "\[Phi]", ",", "\[Mu]1", ",", "\[Mu]"}], 
            "]"}]}], "+", 
          SuperscriptBox["p", "2"]}], "]"}], "]"}], "*", 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["k", "2"], "+", 
          SuperscriptBox["q1", "2"], "+", 
          SuperscriptBox["p", "2"], "-", 
          RowBox[{"2", "*", "k", "*", "q1", "*", "\[Mu]1"}], "+", 
          RowBox[{"2", "*", "k", "*", "p", "*", "\[Mu]"}], "-", 
          RowBox[{"2", "*", "q1", "*", "p", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]1", ",", "\[Phi]", ",", "\[Mu]1", ",", "\[Mu]"}], 
            "]"}]}]}], "]"}], "]"}], "*", 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["p", "2"], "+", 
          SuperscriptBox["q1", "2"], "+", 
          SuperscriptBox["q2", "2"], "-", 
          RowBox[{"2", "*", "p", "*", "q1", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]", ",", "\[Phi]1", ",", "\[Mu]", ",", "\[Mu]1"}], 
            "]"}]}], "-", 
          RowBox[{"2", "*", "p", "*", "q2", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]", ",", "\[Phi]2", ",", "\[Mu]", ",", "\[Mu]2"}], 
            "]"}]}], "+", 
          RowBox[{"2", "*", "q1", "*", "q2", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]1", ",", "\[Phi]2", ",", "\[Mu]1", ",", "\[Mu]2"}],
             "]"}]}]}], "]"}], "]"}]}], ")"}], "/", "64."}], "/", 
    SuperscriptBox["\[Pi]", "3"]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5920663919167843`*^9, 3.592066396532892*^9}, {
   3.592066458555657*^9, 3.592066805894929*^9}, 3.592067410093257*^9, {
   3.592067448640141*^9, 3.592067753453969*^9}, {3.592067798100644*^9, 
   3.5920679169902067`*^9}, {3.592068001218309*^9, 3.592068040062264*^9}, 
   3.592069337368183*^9, {3.59206945968678*^9, 3.592069460103941*^9}, {
   3.5920695056602287`*^9, 3.5920695595299063`*^9}, {3.592069608487071*^9, 
   3.592069751026957*^9}, {3.592150839910246*^9, 3.59215084780333*^9}, {
   3.592153537485997*^9, 3.592153566789893*^9}, {3.592223259730954*^9, 
   3.592223283894908*^9}, {3.5922255267657967`*^9, 3.5922255284008102`*^9}, {
   3.5927612547253933`*^9, 3.592761335605227*^9}, {3.592761494875949*^9, 
   3.592761500818406*^9}, {3.592762052046629*^9, 3.592762121015386*^9}, {
   3.5927622925193157`*^9, 3.592762294594532*^9}, {3.592763655449903*^9, 
   3.5927638177349157`*^9}, {3.592764111914147*^9, 3.592764143295429*^9}, {
   3.5927710299524918`*^9, 3.592771045323942*^9}, {3.592821990166234*^9, 
   3.592821999012002*^9}, {3.59282217043034*^9, 3.592822174090159*^9}, 
   3.59282628851263*^9, {3.592827202419595*^9, 3.5928272029499607`*^9}, {
   3.592827556201213*^9, 3.5928275568993597`*^9}, {3.59282762194914*^9, 
   3.592827650945356*^9}, {3.592850891018014*^9, 3.592850918549065*^9}, {
   3.592851016030954*^9, 3.5928510275058603`*^9}, {3.592851057811228*^9, 
   3.592851118888485*^9}, 3.5928514483024893`*^9, {3.592869530552868*^9, 
   3.592869554335053*^9}, {3.5928695886175547`*^9, 3.592869708033636*^9}, {
   3.592871608495091*^9, 3.5928716798080177`*^9}, {3.593191589501388*^9, 
   3.593191635049176*^9}, {3.593191724405734*^9, 3.593191742744581*^9}, {
   3.593191837638856*^9, 3.5931919100235767`*^9}, {3.593192155329774*^9, 
   3.593192156225503*^9}, {3.593192232513816*^9, 3.593192262598289*^9}, {
   3.593192295004855*^9, 3.593192419986742*^9}, {3.593192490670388*^9, 
   3.593192493915887*^9}, {3.593192534582417*^9, 3.5931925356987762`*^9}, {
   3.59319856535787*^9, 3.593198582837366*^9}, {3.5932197161334267`*^9, 
   3.5932197236399603`*^9}, 3.593276882937324*^9, {3.593387291530457*^9, 
   3.5933872998572598`*^9}, {3.593387330264431*^9, 3.593387340145187*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Mass", "[", "R_", "]"}], ":=", " ", 
   RowBox[{
    FractionBox["4", "3"], "\[Pi]", "*", 
    SuperscriptBox["R", "3"], "*", 
    FractionBox[
     RowBox[{"(", 
      RowBox[{"9.648", " ", 
       SuperscriptBox["10", "11"], "*", 
       RowBox[{
        RowBox[{"(", "1024", ")"}], "^", "3"}]}], ")"}], 
     SuperscriptBox["2400", "3"]]}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.5920695847320747`*^9, {3.5923063497844067`*^9, 3.592306400379266*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"fI21R", "[", "R", "]"}], "]"}], ")"}], "*", 
   RowBox[{"(", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"Sqrt", "[", 
          RowBox[{"fI21R", "[", "RR", "]"}], "]"}], ",", " ", "RR"}], "]"}], "/.",
        " ", 
       RowBox[{"RR", "\[Rule]", "R"}]}], ")"}], 
     RowBox[{"-", "1"}]], ")"}], "*", 
   RowBox[{"(", 
    FractionBox[
     RowBox[{
      RowBox[{"-", "2.23"}], " ", 
      SuperscriptBox["10", "15"]}], 
     RowBox[{"1.43", " ", 
      SuperscriptBox["10", "16"]}]], ")"}]}], "-", "2"}]], "Input",
 CellChangeTimes->{{3.59286996614828*^9, 3.592870006626381*^9}, {
  3.592870053588168*^9, 3.592870110807674*^9}, {3.59287020636482*^9, 
  3.592870253993087*^9}, {3.592870345690892*^9, 3.592870384856512*^9}, {
  3.5928704617709*^9, 3.5928704999916677`*^9}, {3.592870591631545*^9, 
  3.59287060209901*^9}, {3.592870730708084*^9, 3.592870741248386*^9}, {
  3.592870776527237*^9, 3.5928707879290447`*^9}, {3.5928708604416018`*^9, 
  3.592870864113789*^9}, {3.5928712330763073`*^9, 3.592871233604826*^9}, {
  3.592871350263422*^9, 3.592871355717187*^9}, {3.592871438026409*^9, 
  3.592871488891148*^9}, {3.592873820811932*^9, 3.5928738305527363`*^9}, {
  3.5929113941030693`*^9, 3.59291139459459*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "2"}], "-", 
  FractionBox[
   RowBox[{"0.3118881118881119`", " ", 
    RowBox[{"fI21R", "[", "R", "]"}]}], 
   RowBox[{
    SuperscriptBox["fI21R", "\[Prime]",
     MultilineFunction->None], "[", "R", "]"}]]}]], "Output",
 CellChangeTimes->{
  3.592870060909154*^9, {3.592870091465899*^9, 3.592870111066575*^9}, 
   3.592870219902481*^9, 3.5928702549016027`*^9, 3.592870385310302*^9, {
   3.5928704775596867`*^9, 3.592870501387368*^9}, 3.5928706023918962`*^9, 
   3.592870742354405*^9, 3.5928707883379602`*^9, 3.592870864606748*^9, 
   3.592871234028789*^9, {3.592871350917947*^9, 3.592871356300364*^9}, {
   3.5928714675899487`*^9, 3.592871489487945*^9}, 3.5928738315133*^9, 
   3.592914757519402*^9, 3.5931914799122667`*^9, 3.593385029426218*^9, 
   3.5933875462040462`*^9, 3.5934438221226873`*^9, 3.5934439817695217`*^9, 
   3.5934450337142344`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"fI21R", "[", "R", "]"}], "+", 
     RowBox[{
      SuperscriptBox[
       RowBox[{"(", "3000", ")"}], "2"], "*", 
      RowBox[{"fI22R", "[", "R", "]"}]}]}], "]"}], ")"}], "*", 
  RowBox[{"(", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{"fI21R", "[", "RR", "]"}], "+", 
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["3000", "2"], "*", 
            RowBox[{"fI22R", "[", "R", "]"}]}], ")"}]}], "]"}], ",", " ", 
        "RR"}], "]"}], "/.", " ", 
      RowBox[{"RR", "\[Rule]", "R"}]}], ")"}], 
    RowBox[{"-", "1"}]], ")"}], "*", 
  RowBox[{"(", 
   FractionBox[
    RowBox[{
     RowBox[{"-", "1.3657"}], " ", 
     SuperscriptBox["10", "15"]}], 
    RowBox[{"7.999", " ", 
     SuperscriptBox["10", "15"]}]], ")"}]}]], "Input",
 CellChangeTimes->{{3.5928712049493628`*^9, 3.592871216311335*^9}, 
   3.592871745969755*^9, {3.592873850289229*^9, 3.592873900119732*^9}, {
   3.592910755048575*^9, 3.592910766733501*^9}, {3.592910874385602*^9, 
   3.59291089182335*^9}, {3.5929112903922*^9, 3.592911343439502*^9}, {
   3.592911398498057*^9, 3.592911462320829*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"0.34146768346043255`", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"fI21R", "[", "R", "]"}], "+", 
      RowBox[{"9000000", " ", 
       RowBox[{"fI22R", "[", "R", "]"}]}]}], ")"}]}], 
   RowBox[{
    SuperscriptBox["fI21R", "\[Prime]",
     MultilineFunction->None], "[", "R", "]"}]]}]], "Output",
 CellChangeTimes->{{3.592911414030724*^9, 3.592911463311989*^9}, 
   3.5929147576752253`*^9, 3.593191479989336*^9, 3.593385029492236*^9, 
   3.593387546269198*^9, 3.59344382279175*^9, 3.593443982159513*^9, 
   3.593445034317615*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"fI21R", "[", "R_", "]"}], ":=", 
   RowBox[{"fI21", "[", 
    RowBox[{"Mass", "[", "R", "]"}], "]"}]}], ";", 
  RowBox[{
   RowBox[{"fI22R", "[", "R_", "]"}], ":=", 
   RowBox[{"fI22", "[", 
    RowBox[{"Mass", "[", "R", "]"}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.592871418756938*^9, 3.592871433896693*^9}, {
  3.592911322330215*^9, 3.5929113317961206`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fI21R", "[", "RR", "]"}]], "Input",
 CellChangeTimes->{{3.59291135770126*^9, 3.592911390196335*^9}, {
  3.592918240010529*^9, 3.592918248269751*^9}}],

Cell[BoxData[
 RowBox[{"fI21", "[", 
  RowBox[{"3.139005299315886`*^11", " ", 
   SuperscriptBox["RR", "3"]}], "]"}]], "Output",
 CellChangeTimes->{{3.592918244305825*^9, 3.592918248542329*^9}, 
   3.5931914800984907`*^9, 3.593385029575032*^9, 3.5933875463162117`*^9, 
   3.593443824338129*^9, 3.593443982683691*^9, 3.5934450345874*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"SigmaRSq", "[", 
   RowBox[{"q_", ",", "ftNL_", ",", "R_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"fI21R", "[", "R", "]"}], "+", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"q", "*", "ftNL"}], ")"}], "2"], "*", 
    RowBox[{"fI22R", "[", "R", "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.5929183044555893`*^9, 3.5929183156174707`*^9}, {
  3.592918449308448*^9, 3.592918472481195*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SigmaRSq", "[", 
  RowBox[{"q", ",", "ftNL", ",", "RR"}], "]"}]], "Input",
 CellChangeTimes->{{3.5929184803356457`*^9, 3.5929184887579193`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"fI21", "[", 
   RowBox[{"3.139005299315886`*^11", " ", 
    SuperscriptBox["RR", "3"]}], "]"}], "+", 
  RowBox[{
   SuperscriptBox["ftNL", "2"], " ", 
   SuperscriptBox["q", "2"], " ", 
   RowBox[{"fI22", "[", 
    RowBox[{"3.139005299315886`*^11", " ", 
     SuperscriptBox["RR", "3"]}], "]"}]}]}]], "Output",
 CellChangeTimes->{3.592918488981443*^9, 3.593191480183752*^9, 
  3.593385029604027*^9, 3.593387546364791*^9, 3.5934438268820543`*^9, 
  3.5934439828296137`*^9, 3.593445035208312*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"D", "[", 
  RowBox[{
   RowBox[{"SigmaRSq", "[", 
    RowBox[{"q", ",", " ", "ftNL", ",", " ", "RR"}], "]"}], ",", " ", "RR"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.592918638485166*^9, 3.592918647064118*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"9.417015897947659`*^11", " ", 
   SuperscriptBox["RR", "2"], " ", 
   RowBox[{
    SuperscriptBox["fI21", "\[Prime]",
     MultilineFunction->None], "[", 
    RowBox[{"3.139005299315886`*^11", " ", 
     SuperscriptBox["RR", "3"]}], "]"}]}], "+", 
  RowBox[{"9.417015897947659`*^11", " ", 
   SuperscriptBox["ftNL", "2"], " ", 
   SuperscriptBox["q", "2"], " ", 
   SuperscriptBox["RR", "2"], " ", 
   RowBox[{
    SuperscriptBox["fI22", "\[Prime]",
     MultilineFunction->None], "[", 
    RowBox[{"3.139005299315886`*^11", " ", 
     SuperscriptBox["RR", "3"]}], "]"}]}]}]], "Output",
 CellChangeTimes->{3.5929186473331833`*^9, 3.593191480256428*^9, 
  3.593385029628903*^9, 3.593387546388138*^9, 3.593443827293437*^9, 
  3.593443982851913*^9, 3.593445035378461*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.593275904263381*^9, 3.593276121372016*^9}, {
   3.5932761547888613`*^9, 3.593276323502252*^9}, 3.593277003625215*^9, {
   3.5932770414094048`*^9, 3.5932770736175423`*^9}, {3.5932772571596622`*^9, 
   3.593277288587118*^9}, 3.5933873261290913`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Ik22Test", "[", 
  RowBox[{
  "0.03", ",", "R", ",", "0.67217", ",", " ", "0.5", ",", " ", "1.5708", ",", 
   " ", "2.011", ",", " ", 
   RowBox[{"-", "0.5"}], ",", " ", "4.71239", ",", " ", "0.67217", ",", " ", 
   "0.5", ",", " ", "4.71239"}], "]"}]], "Input",
 CellChangeTimes->{{3.59327633024555*^9, 3.5932763631840687`*^9}, {
  3.5932767760675907`*^9, 3.593276860708712*^9}}],

Cell[BoxData[
 RowBox[{"Ik22Test", "[", 
  RowBox[{
  "0.03`", ",", "R", ",", "0.67217`", ",", "0.5`", ",", "1.5708`", ",", 
   "2.011`", ",", 
   RowBox[{"-", "0.5`"}], ",", "4.71239`", ",", "0.67217`", ",", "0.5`", ",", 
   "4.71239`"}], "]"}]], "Output",
 CellChangeTimes->{{3.593276364760063*^9, 3.593276366587741*^9}, 
   3.5932768610411777`*^9, 3.5932768915871477`*^9, 3.593277007400895*^9, 
   3.59338502967032*^9, 3.593387546418049*^9, 3.593443827903203*^9, 
   3.5934439830094967`*^9, 3.5934450355553303`*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{609, 721},
WindowMargins->{{Automatic, 20}, {55, Automatic}},
Magnification:>0.75 Inherited,
FrontEndVersion->"9.0 for Linux x86 (64-bit) (February 7, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 298, 4, 27, "Subsubsection"],
Cell[880, 28, 22173, 654, 151, "Input",
 InitializationCell->True],
Cell[CellGroupData[{
Cell[23078, 686, 197, 4, 2, "Input"],
Cell[23278, 692, 321, 4, 70, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23636, 701, 494, 12, 3, "Input"],
Cell[24133, 715, 26518, 551, 70, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[50688, 1271, 216, 5, 2, "Input"],
Cell[50907, 1278, 425, 11, 70, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[51381, 1295, 327, 4, 29, "Subsection"],
Cell[51711, 1301, 19160, 518, 1046, "Input"],
Cell[70874, 1821, 504, 15, 44, "Input"],
Cell[CellGroupData[{
Cell[71403, 1840, 1378, 33, 65, "Input"],
Cell[72784, 1875, 892, 18, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[73713, 1898, 1295, 37, 85, "Input"],
Cell[75011, 1937, 595, 15, 42, "Output"]
}, Open  ]],
Cell[75621, 1955, 417, 11, 24, "Input"],
Cell[CellGroupData[{
Cell[76063, 1970, 173, 3, 24, "Input"],
Cell[76239, 1975, 336, 6, 28, "Output"]
}, Open  ]],
Cell[76590, 1984, 440, 12, 26, "Input"],
Cell[CellGroupData[{
Cell[77055, 2000, 168, 3, 24, "Input"],
Cell[77226, 2005, 529, 13, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[77792, 2023, 237, 6, 24, "Input"],
Cell[78032, 2031, 803, 20, 47, "Output"]
}, Open  ]],
Cell[78850, 2054, 296, 4, 24, "Input"],
Cell[CellGroupData[{
Cell[79171, 2062, 404, 8, 42, "Input"],
Cell[79578, 2072, 518, 10, 42, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
